Quantum bistability and spin current shot noise of a single quantum dot coupled to 

an optical microcavity 
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Here we explore spin dependent quantum transport through a single quantum dot coupled to an 
optical microcavity. The spin current is generated by electron tunneling between a single doped 
reservoir and the dot combined with intradot spin flip transitions induced by a quantized cavity 
mode. In the limit of strong Coulomb blockade, this model is analogous to the Jaynes-Cummings 
model in quantum optics and generates a pure spin current in the absence of any charge current. 
Earlier research has shown that in the classical limit where a large number of such dots interact 
with the cavity field, the spin current exhibits bistability as a function of the laser amplitude that 
drives the cavity. We show that in the limit of a single quantum dot this bistability continues to 
be present in the intracavity photon statistics. Signatures of the bistable photon statistics manifest 
themselves in the frequency dependent shot noise of the spin current despite the fact that the 
quantum mechanical average spin current no longer exhibits bistability. Besides having significance 
for future quantum dot based optoelectronic devices, our results shed light on the relation between 
bistability, which is traditionally viewed as a classical effect, and quantum mechanics. 

PACS numbers: 42.50.Pq,73.63.Kv,78.67.Hc 



I. INTRODUCTION 

Bistability is a phenomenon that readily occurs in clas- 
sical systems that possess a nonlinear response to some 
input signal. In a bistable system the output function, 
F{I), can exhibit two stable states for a certain range 
of the input / such that when / is varied F{I) follows a 
hysteresis loop. One of the most familiar examples is the 
hysteresis curve in the magnetization of a ferromagnetic 
material in the presence of an external magnetic field. 
In the context of electronics, digital flip-flop circuits and 
Schmitt triggers are common examples of bistable cir- 
cuits. In nonlinear optics, optical bistability (OB) occurs 
in the input-output function of an optical resonator that 
contains a nonlinear dielectric and is driven by a laser—. 
OB has a number of applications in optical communi- 
cations and computing because it can be used to build 
all optical switches, logic gates, and optically bistable 
memory devices^i^i^i^i^ but is also interesting for basic 
studies of phase transitions between stationary but non- 
equilibrium states^. 

Here we explore a model first proposed by two of ns^'^ 
that unifies research in nonlinear quantum optics with 
spintronics. In the present work, we use that model to 
explore how bistability manifests itself in the quantum 
world. Spintronics has emerged as a field in which the 
spin degrees of freedom of charge carriers in solid state 
devices are exploited for the purpose of information pro- 
cessing. Manipulation of the spin degrees of freedom 
rather than the charge has the advantage of longer coher- 
ence and relaxation times since the spin is more weakly 
coupled to its environment^^. For the same reason, ma- 
nipulation of the spin of an electron is much harder than 
the charge and therefore has resulted in significant effort 
to come up with proposals for necessary spin devices in- 
cluding spin batteries, spin filters, spin transistors, etc... 



Much of this work has focused on ways to generate pure 
spin currents, I^ = s(/| — Ii)i which are the result of an 
equal number of spin up (t) and spin down (|) charge car- 
riers moving in the opposite direction so that the charge 
current, Ic = q{I] + -^x), is zero. Here, !„ are the spin 
polarized particle currents, s = Ti/2 the spin of the parti- 
cle, and q = e the charge. There currently exist numerous 
theoretical and experimental concepts for generating spin 
currents in semiconductor nanostructures including spin- 
orbit (SO) interactionsiiii^, optical absorption^'^ and Ra- 
man scatteringi^, as well as various types of quantum 

Electron spin resonance (ESR) between Zeeman states 
in a quantum dot connected to leads is one of 
the proposed models for the generation of pure spin 
currents^ijS . According to this model, spin flips are the 
result of a transverse magnetic field that cause the spin 
direction of outgoing electrons to be opposite to those 
entering the dot. Our model^'^ extends this idea to spin 
flips induced by Raman transitions inside of an optical 
microcavity. One laser involved in the Raman transition 
is a strong undepleted pump while the other is a mode of 
the cavity. Inside of the cavity, both the feedback effect 
resulting from light "bouncing" back and forth numerous 
times in the cavity and quantum fluctuations can have a 
dramatic influence on the characteristics of the spin cur- 
rent. In our previous work^, we considered the classical 
limit of a large number of dots, N ^ 1, interacting with 
the cavity mode such that quantum noise is negligible. 
When the cavity is driven by a laser, the system exhibits 
absorptive OB in the amplitude of the cavity field. Be- 
cause the spin current is a function of the cavity field 
amplitude, the spin current also exhibits bistability as 
function of the amplitude of the driving laser which sur- 
vives even in the presence of significant variations in the 
dot sizes and coupling to the cavity field. 



However, this bistability is a purely classical effect 
since a large number of dots collectively interact with 
the cavity mode like a single classical absorber. This 
begs the question of what happens if we consider only 
a single quantum dot coupled to the cavity where quan- 
tum fluctuations will be so large as to imply that the 
two 'stable' outputs lose their stability. Earlier theo- 
retical work in quantum optics that explored the limit 
of 'bistability' for a single atom coupled to cavity mode 
found that the steady state phase space distribution of 
the cavity field had a bimodal structure indicative of two 
'stationary' value a^'^i^'^i^^'^^ . These states are not how- 
ever stable since quantum noise forces stochastic jumps 
between the two values^. Here we show that, while the 
average spin current for a single dot in a driven cavity 
does not exhibit bistability, the frequency dependent spin 
current shot noise does exhibit signatures of the two sta- 
tionary cavity states since in this system the shot noise 
spectrum reflects the probability distribution for cavity 
photon states. In contrast to the mentioned work from 
quantum optic3^iSii^i2£, which relied on quantum tra- 
jectory Monte Carlo simulations, we utilize a standard 
master equation to calculate the shot noise indicating 
that evidence of the quantum limit of bistability can be 
gleaned using more pedestrian techniques. 

In Section II, we briefly review our model and intro- 
duce our mathematical formulation of the shot noise in 
terms of the dot+cavity master equation. In Section HI, 
we numerically study both the average spin current and 
the associated shot noise. In Section IV, we present our 
conclusions. 



II. MODEL 
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FIG. 1: (a) Raman transition between the dot Zeeman states, 
I t) i)i via an intermediate trion state, \+t), induced by a laser 
with frequency lul and a cavity mode with frequency ujcav- 
Both optical fields are detuned from the trion energy by the 
amount Ar implying that | + f ) is a virtual state. The spin 
eigenstates along the direction of the magnetic field are super- 
positions of spin eigenstates in the growth direction, cl^^|0). 
(b) Schematic of a single quantum dot indicating Zeeman en- 
ergy levels in the dot and allowed tunneling between the lead 
and dot. Also shown to the right is a hypothetical configura- 
tion of a dot in a micropillar cavity showing the direction of 
the magnetic field, pump laser, and cavity decay. 



We consider a self-assembled quantum dot embedded 
in a high-Q microcavity, as depicted in Fig. 1. We are 
interested in simultaneous coupling of a dot to a cav- 
ity mode and electrical transport through the dot due 
to tunneling from a doped reservoir. A number of ex- 
periments have already measured the conductance and 
shot noise through individual self-assembled quantum 
dots^>^>^>^ as well spectroscopy of exciton and charged 
exciton states in quantum dots with controllable charg- 
ing from a doped lead^ii^i^^. Other experiments have 
demonstrated strong coupling of individual dots to a sin- 
gle optical microcavity mode^ii^. Recently several of 
these directions have come together in the experiment 
by Strauf et al^ showing a high efficiency single pho- 
ton quantum dot source. The experiment demonstrated 
electrical gate controlled charging of dot, which was em- 
bedded in a high-Q optical microcavity, from an n doped 
layer. Several other experiments have followed demon- 
strating electrically driven quantum dots embedded in 
high-Q micropillar cavities that behave as single photon 
sources^^i^. 

We assume that a single electron reservoir at chem- 
ical potential, /^, is coupled to the dot via tunneling. 



Only a single empty orbital energy level, e, of the dot 
lies close to /i. The Zeeman splitting between the two 
electron spin states is A = e^ — £| = Qx^ibB where B 
is a static magnetic field along the x-axis that is per- 
pendicular to the growth direction (z). /zb is the Bohr 
magneton and g^ is the electronic g-factor along the di- 
rection of the magnetic field. The energy levels satisfy 
e^ — e — A/2 < /i < £| = e + A/2 so that only spin up 
electrons can tunnel into the dot and only spin down elec- 
trons can tunnel out. In the limit of very large Coulomb 
blockade energy, which we consider here, only a single 
electron from the reservoir can occupy the dot. The Zee- 
man states along the direction of the B field are super- 
positions of spin eigenstates along the growth direction, 
I t, i) = (1/^2) (cj JO) ± c\^ 10)] , where c^ is an electron 
creation operator. 

Raman transitions between the dot Zeeman states, 
I tii)? via an intermediate trion state, | -f t), are in- 
duced by a a^ polarized laser with frequency uj^ and a 
linearly polarized cavity mode with frequency tOcav Sev- 
eral experiments have already demonstrated the use of 
Raman scattering via an intermediate trion state to ma- 



nipulate electron spin states in quantum dot o^^i'^^'^^i'^^ 
and theoretically such processes have been studied inside 
of optical microcavities for use as a quantum computer—. 
The a~^ pump creates a +3/2 heavy hole and an elec- 
tron with spin down along the z direction according to 
the Hamiltonian, Hpump — {fini/2) exp{—iujit)c] h,^/2j 
which couples to the component of the dot Zeeman states 
with spin up along z yielding a trion state with an elec- 
tron singlet. The a^ component of the cavity field along 
with the pump leads to Raman transitions via the in- 
termediate I -I- 1) state that flips the electron spin while 
the a^ component gives rise to additional energy shifts 
due to the AC stark effect. When the two fields are far 
detuned by an amount A^ from the creation energy for 
the I + 1) state, the intermediate trion state can be adi- 
abatically eliminated to give Hi = ^(aclcfe*'^'* — /i.e.) 
where g = gcav^i/'^^R and a is the photon annihilation 
operator for the cavity mode^i^. We have absorbed all 
energy shifts of the states \a) due to the AC stark ef- 
fect into a renormalization of the energy levels £„■ Non- 
resonant terms atcjcie^*"'* + h.c. can be neglected pro- 
vided |A - {uJcav - W;)| < |A -I- {uOcav - ^l)\- 

As one can see in Fig. 1, if an electron enters the dot 
in the spin "f state, a photon must be absorbed from 
the cavity mode and emitted into the pump in order 
to generate a spin current. It is therefore necessary to 
drive the cavity field. We assume that the cavity is 
driven by a classical source oscillating at frequency Wp, 
Hp — ihe{eyip{—iijjpt)a^ — h.c), corresponding to coher- 
ent coupling between a laser and the cavity modei>S. 

The Hamiltonian in a frame rotating at the frequency 
LjpisH' = H^+H'p + H'j, 



K 



h{w,av - iOp)A^A + e{C]C^ + CJCx) 
+ {\ + u:i~Up){C\C^-C\C^)/2 (I) 
H'j + H'p = ihg{AClC^ - h.c.) + ih€{A^ - h.c.) (2) 

Here, we have defined operators in a rotating frame 
a = Ae^'^"*, C| = (7| exp(— i(wi — u}p)t/2), and C| = 
Ci exp(i(a;/ — LUp)t/2). In this work we assume that the 
resonance conditions, Ucav — ^p and A = cjp — wj, are al- 
ways satisfied, so that the final Hamiltonian of the system 
is H' = e{C\C^^+ClCi) + ihg{AClC^ - h.c.) + ih€{A^ - 
h.c.) 

The dynamics of the system can be described in terms 
of the density operator, p, for the cavity plus dot. The 
master equation for p is given by, 

p = -i[H\p\lh-T,av{A^ Ap^2ApA^ + pA^ A)/2 + p\i,ad 

(3) 
The first term describes coherent dynamics of the cou- 
pled QD-cavity system, the second term represents the 
cavity decay-^'^'^, and the third term describes QD-lead 
coupling. The lead-dot coupling is most easily expressed 
in terms of the matrix elements of the density opera- 
tor, pj^^ = {n, a\p\a' , m) where |cr, n) represents a state 



with n photons in the cavity and cr = 0, t, i correspond- 
ing to no electrons, one spin up, or one spin down elec- 
tron, respectively. The specific form of the master equa- 
tions for the lead coupling are^i^^ 



. (n,77i) I 
Po,0 \lead 

. (n,m) I 
Pt t \lead 

. (n,m) I 
P\ I \lead 

. (n,m) I 
Pt I \lead 



p(_) (n,m) _-p(+) ( 
^ Pi,i ^ ^" 



Po,o 



p(+) («,m) 



(4) 
(5) 
(6) 
(7) 



Here, r*-^-' is the rate at which spin down electrons tun- 
nel out of the dot into lead and F^^^ is the rate at which 
spin up electrons tunnel into the dot. We assume that 
the tunnelling between the lead and the dot is spin in- 
dependent, F'^^' = F'^' = F. We can rewrite Eq. [3] in 
matrix form, 



dp/dt = Mp 



(8) 



T (n:m) 

where p , 



p is the density matrix in vector form. 

The steady state solution, p, is given by the eigenvector 
of M with zero eigenvalue. Conservation of probability 
insures that M has a zero eigenvalue^. 

The spin current operator is defined as, Is = s(/| —li) 
with the stationary currents given by (/|) = Fpo,o and 

(A) = -rPi,i- Here /5^^a' = E„ pi"^"' is the reduced 
density matrix of the dot after tracing over the cavity 
field and s = h/2. We note that the spin current can be 
easily interpreted as the rate at which spin up electrons 
tunnel into the empty dot, Fpo.o, plus the rate at which 
spin down electrons leave the dot, — Fpj^j^. The average 
spin current can be expressed in terms of expectation 
values of the cavity field using Eq. [31 



{Is)^2s{2eRe[{A)]^T,av{A^A)) 



(9) 



One sees that the spin current is also the difference be- 
tween the rate at which photons are coherently injected 
into the cavity by the driving laser, 2eRe[{A)], and the 
rate at which photons decay from the cavity, Fcau(^^^)- 
Conservation of energy, implies that this difference must 
be absorbed by a spin flip of the electron in the dot. 

The noise power spectrum for the current can be ex- 
pressed as the Fourier transform of the current-current 
correlation function, 



Sa,a'iuj) 



dte'"*[(j,(t)j,Ko)) -(/.)(/.')]• (10) 



The spin current shot noise, S*^^^ = 

2 j^^dtexp{iujt)[{isit)is{0)) ~ {Is) (Is)] can be written 
in terms of the shot noise spectrum for the spin resolved 
currents as S'^^^ = s^(5'|,-f + S'j,,^ - S^^i - Si^i). It is 
well known that for currents comprised of uncorrelated 
particles, the noise power spectrum is Poissonian, 
S{uj) = 2q{I), where q is the quantity transported by 



each particle in the current J^s^ 5 = e in the case of 
standard charge currents while in our case the trans- 
ported quantity is spin q = s. It is often convenient to 
measure the shot noise relative to the Poissonian noise 
by defining the Fano factor, 



FH 



5(^)(lj) 



2s/, 



(11) 



where Is is the average spin current. F{uj) > 1 represents 
super-Poissonian noise while F{uj) < 1 represents sub- 
Poissonian noise. 

Here we adopt the numerical method for evaluating 
Eq. [To] developed in Ref.— for use with master equations 
of the form Eq. [5] Briefly stated, the spectral decompo- 
sition of the matrix M is given by M = ^^^ XP\ where A 
is an eigenvalue of M and P\ is the projection operator 
associated with that eigenvalue. This form of M can be 
used to evaluate the time evolution of the current opera- 
tors, la, and in the end yields the following form for the 
spin current shot noise. 

Si^\.) ^ 2sls + 2 V M^P^if + Tr[LPJp]\ 

(12) 

and the first term, the Poissonian contribution, is cal- 
culated from Is = Tr[Isp\- Here we note that the 
projection operators can be calculated in terms of the 
left and right eigenvectors of M, P\ = v\ [wx)^ where 
{wx) M = A(wa) define the left eigenvectors while 
Mv\ = \v\ define the right eigenvectors. They satisfy 
the orthonormality relation {wx^y vxm — ^n,m- We note 
that this is a different formulation of the projection op- 
erators than appears Ref4^ where the Px = SEnS~^, 
S being the matrix whose columns are the right eigen- 
vectors of M and En is a square matrix that has zero 
entries everywhere except the (n, n) element, which is 1. 
However, it is easy to show that these forms are mathe- 
matically equivalent. 



III. RESULTS 

We first discuss the behaviour of the intracavity field 
as a function of the driving field amplitude, e. The in- 
tracavity field can be readily visualized in term of the 
Q-distributionii^ for the cavity mode in the steady state 
as shown in Fig. 2. The Q-distribution is defined as 



Q(a) 



X! (a,cr|p|a,cr)/7r 
■=o,T,i 



where \a) is a coherent state A\a) = a\a). It represents a 
pseudo-quantum mechanical phase space distribution for 
bosonic quantum fields where _Re[a] and /m[Q;], which 
represent the quadrature components of the field, can be 
interpreted as the position and momentum, respectively. 



of a fictitious particle. We note that there exist a number 
of different pseudo-phase space distributions for bosonic 
fields whose utility depends on the particular problem^. 
We chose the Q-distribution because it is both positive 
semi-definite and can also be interpreted as a probabil- 
ity distribution, namely the probability of measuring the 
field in the coherent state \a). This therefore allows qual- 
itative comparisons to classical phase space probability 
distributions. 

In Fig. 2(a), which corresponds to weak driving, there 
is only single peak around a « 0. This represents a cav- 
ity that is over damped such that all energy injected into 
the cavity is absorbed by the dot. For larger driving, as 
in Fig. 2(b) and (c), there are two peaks, one at a ~ 
and another at i?e[a] > and /to [a] = 0. This repre- 
sents the bistable situation where the cavity field has two 
most probable states. By contrast, in Fig. 2(d) one can 
see that the peak around a — has completely disap- 
peared and only a peak with i?e[a] > remains when 
the driving is further increased. This peak corresponds 
to the case where the cavity driving is so strong that the 
dot transition is saturated. For a saturated transition, 
Po.o — Pt,T ~ Pi.i ~ 1/3 such that the current obtains 
the maximum value, (/g) = s(rpo,o + ^Pi.l) = 2sr/3. 
Based on Eq. [51 the approximate location of this second 
peak in the Q-distribution is then 



(2e/r,,,) -I- V(2e/rca.)2 - 4r/3r, 



(13) 



where we note that the last term due to the lead, 
AT/STcav, reduces the cavity field amplitude below the 
value of an empty cavity (i.e. no absorber in the cav- 
ity), 2e/rcav Fig- 3 shows the peak values of the Q- 
distribution as a function of e where one can see that a 
classic hysteresis loop emerges. This can be compared to 
the semiclassical solution for the cavity amplitude that 
ignores quantum fluctuations. 



|e| -^cav\ 



gM^ 



6<?2|a|2-Kr2/2 



(14) 



Equation [T4| is obtained from the equations of motion for 
the expectation values of the cavity and dot operators 
by factorizing the expectation values of products of op- 
erators such as (A^CJCi) -^ {A^){C^^Ci) (Note that one 
recovers Eq. [T^jfrom Eq. [Hjin the limit that g\a\ ^ F.). 
One can see in Fig. 3 that in the quantum case, the range 
of e values where bistability is present has been shifted 
to higher values due to quantum fluctuations. 

By contrast, in Fig. 4, we present the quantum me- 
chanical average spin current for a single dot as a func- 
tion of the driving amplitude. As one can see it is a single 
valued quantity that shows no sign of the 'switch back' 
behaviour characteristic of bistability that is seen in the 
inset, which is the semiclassical spin current calculated 
using Eq. [Hjand Is = 2s{2\a\\e\-Tcav\a\'^)- In fact, the 
current is qualitatively the same as that calculated for 
spin flips in the case of ESR using a classical magnetic 



o^Q (a 




0.7 

0.6 

0.5 

.0.4 



- 0.3 



0.2 



0.1 







/ 0.8 

/ ■ 0-8 

/ I"* 
/ J> 0.2 




- 




- 


\ 


/ % 0.2 0.4 0.6 0.8 

_^ — III 





0.2 



Irti(u) 



0.4 06 

8 



08 



FIG. 2: (Color Online) Q-distribution vs. _Re[ci] and /mfa] 
for Fcav = 0.2r, g = 2r and in clockwise order (a) e — 0.43r 
(b)e = 0.53r (c) e = 0.55r, (d) e = 0.7r 
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FIG. 3: (Color Online) Peak values of the Q-distribution as a 
function of e (in units of F) for Vcav = 0.2r, g — lAT (circles) 
and Vcav = 0.2r, and g = 2F (diamonds). For comparison, 
the semiclassical solution Eq. [TJ] exhibiting bistability is also 
shown for the same parameters (green and black solid lines). 
One can see that the Q-distribution for a single dot qualita- 
tively follows the semiclassical solution although the range of e 
where 'bistability' occurs is reduced by quantum fluctuations. 



field^-^ . This is not surprising since one can see from Eq. 
[5] that the spin current is the quantum mechanical expec- 
tation value of the cavity field and despite the bimodal 
distribution of Q{a), the spin current is averaged over 
both values, {Is) ~ Pils,i + P2ls,2 where Pj are the to- 
tal probabilities corresponding to each of the two peaks 
in Q{a) and Isj — 2s{2eRe[aj] — TcavCt*aj) where aj 
are the locations of the two peaks. 

This begs the question, how does the bistable struc- 



FIG. 4: (Color Online) Average value of the spin current, 
Is (in units of sF) as a function of e (in units of F) for 
Fcati = 0.2F, g = 1.4F. Inset shows the semiclassical spin 
current obtained from Eq. I14l and Is = 2s(2|a||e| — Fcai>|a| ). 
One can see that when quantum fluctuations are included, all 
indications of bistability are destroyed in the average current. 



ture of the intracavity field manifest itself in quantum 
mechanical observables? Previous work on the quan- 
tum limit of bistability for single atom cavity QED fo- 
cused on the quantum dynamics using 'quantum tra- 
jectories' Monte Carlo simulations approach based on 
stochastic Schrodinger equations^^ and stochastic mas- 
ter equation o^^'^^ , which showed that the cavity field and 
photocurrent from the cavity undergo stochastic jumps 
between the two states given by the peaks in the Q- 
distribution. In these systems, the average time between 
switching events was proportional to the spontaneous 
emission lifetime since it was the 'wave function collapse' 
due to spontaneous emission of the atom that drove the 
system between the two states^. 

Equations [4][7] have a similar form to that of the mas- 
ter equation for atomic decay. Therefore we can draw an 
analogy with earlier work and argue that 'wave function 
collapse' resulting from electron tunneling events into 
and out of the dot will induce jumps between the two 
stable quantum states of the cavity field. Since the time 
scale that determines transport through the dot is deter- 
mined by the time needed for a spin flip, which is the Rabi 
frequency g\/n, different cavity field states will result in 
different time intervals between successive electrons be- 
ing 'emitted' by the dot. One would therefore expect that 
the Rabi frequencies associated with the two stable field 
states would manifest themselves in the current-current 
correlations, {I{t + T)I{t)). 

Fig. 5 and Fig. 6 show F{w) and P{n), the prob- 
ability distribution for the cavity photons, for different 
values of e. For the sake of comparison. Fig. 7 shows the 
Fano factor for the case of electron spin resonance (ESR) , 
Fesr(w), with a classical field of Rabi frequency R that 
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FIG. 5; (Color Online) (a) Fano factor F{lu) for Fcav = 0.4r, 
g — 2r and different values of e (in units of F). (b) Steady 
state probability distribution for photons in the cavity for the 
same parameters as in (a). Note that ui is measured in units 

of r. 



FIG. 6; (Color Online) (a) Fano factor F{lu) for Fcav = 0.2r, 
g = 1.4r and different values of e (in units of F). (b) Steady 
state probability distribution for photons in the cavity for the 
same parameters as in (a). Note that ui is measured in units 

of r. 



flips the spins of the electrons^-^. The classical field ESR 
Hamiltonian can be obtained by replacing A and A^ with 
a c-number (A — > a) in H' with R = ga. The similarity 
between Eqs. |4][7] and that of spontaneous emission al- 
lows us to define the critical dot number Nq = 2TTcav/g^, 
which represents in the semiclassical theory the minimum 
number of dots necessary for bistability to be present^iS^, 
as well as the critical photon number ric — r^/4(7^, which 
defines the number of photons necessary to significantly 
modify the dot response^^. Classical bistability is pre- 
dicted to occur in the limit Uc —> oo and therefore larger 
values of Uc should produce more pronounced 'bistabil- 
ity' in the single dot/atom case^i. For both Figs. 5 and 
6 No f^ 0.2 while for Fig. 5 ric = 0.063 and for Fig. 6 
Uc = 0.13. This behavior with ric is confirmed in the 
figures where the bimodality of P{n) is more visible and 
present for a larger range of values in Fig. 6 as compared 
to Fig. 5. 

In these figures, we can see that for small e, below 
the threshold for the onset of bistability, F{u!) is super- 
Poissonian for low frequencies and Poissonian at high fre- 
quencies, which is similar to the case of ESR for small 
R where i^esr(O) -^ 3/2 for R ^ 0. For small e, the 
cavity is overdamped and only the vacuum state has sig- 
nificant probability, P(0), and therefore transitions are 
primarily driven by fluctuations above the vacuum state. 



In the opposite extreme with stronger e in the bista- 
bility region, which is most clearly seen in Fig. 6 for 
e/F = 0.42, 0.45, 0.48, F{uj) remains super-Poissonian 
at zero frequency while at a; « 2g\a2\ a broad sub- 
Poissonian dip develops whose overall width is deter- 
mined by the width of P{n) around the second maximum 
at 71,2 = |a2p- This behavior is a mixture of the ESR 
system for small and large R since as already mentioned, 
Fesr is super-Poissonian at low frequencies R <S^ T. By 
contrast, the ESR system exhibits a sub-Poissonian dip 
at 2R for R > T while being nearly Poissonian at zero 
frequency (Fesr(O) -^ 19/18 for R -^ cx)). For even larger 
e such as e/F = 1 in Fig. 5 or e/F — 0.6 in Fig. 6, which 
place the system outside of the bistable regime, one can 
see that the broad sub-Poissonian dip around lu « 2g\a2\ 
persists but that i^(0) is no longer super-Poissonian but 
rather has become sub-Poissonian. Therefore we can con- 
clude that the super-Poissonian behavior of F(0) is at- 
tributable to the maximum in P{n) at n = while the 
sub-Poissonian dip is attributable to the the maximum 
in P{n) around |a2p- 



IV. CONCLUSIONS 
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FIG. 7: (Color Online) Fano factor Fesr('^) for electron spin 
resonance with a classical field of Rabi frequency R that flips 
the spins of the electrons. The zero frequency Fano factor 
is given analytically by Fe,.(0) = ^"^"a+ga+gflt'l''^ ■ The fre- 
quency uj is in units of F and for (a) from top to bottom are 
plotted R/V = 0, R/V = 0.1, R/V = 0.2, R/T = 0.3, and 
R/T — 0.4 while for (b) in order of the minima going from 
left to right R/T = 1, R/F = 2, R/T = 3, R/T = 4, R/T = 5, 
and R/T = 6. As in the previous graphs, lj is in units of F. 



Here we have analyzed the spin current and shot noise 
from a single quantum dot embedded inside of a driven 
optical microcavity. We have shown that as a result of 
the cavity field induced spin flips, the quantum bistability 
present in the cavity field Q-distribution manifests itself 
also in the spin current shot noise from the dot. These re- 
sults indicate that despite the large quantum fluctuations 
that wipe out all trace of the bistability in the average 
current, the shot noise reveals the underlying bimodal 
distribution of the cavity field. This works shows that 
there is no need to make recourse to more complicated 
methods such as stochastic wave function methods in or- 
der to detect bistability in the presence of large quantum 
fluctuations. 
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